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We demonstrate that the RHIC data for hadron production in d-Au collisions for all available rapidities are 
compatible with geometric scaling. In order to establish the presence of scaling violations expected from small- a; 
evolution a much larger range in transverse momentum and rapidity needs to be probed. We show that the fall-off 
of the transverse momentum distribution of produced hadrons at LHC is a sensitive probe of small- a; evolution. 



It is well known that the small-x DIS data show 
geometric scaling [1 . This means that the cross 
section depends on Q^/Q'^{x) only, instead of 
and X independently, where Qs{x) is known as the 
saturation scale. Since geometric scaling arises as 
a feature of saturation from nonlinear evolution 
equations, such as the BK equation [51, its occur- 
rence is often seen as an indication of saturation. 
Here, we investigate this issue by studying the 
scaling properties of hadron collision data from 
RHIC at similarly small x. 

The small- a; inclusive HERA data were shown 
to be well described by the phenomenological 
Golec-Biernat WiisthofT (GBW) model for the 
dipole cross section 3 . In the GBW model, 
the cross section is given by cr = coNQ-Q^ff^rt^x), 
where CTq — 23 mb and the scattering amplitude 
A^GBW is given by 



GBW 



(rt,x) = 1 - exp 



'M) , (1) 



where rt is the dipole size. This amplitude de- 
pends on rlQl{x) only, leading to a Q^/Ql{x) de- 
pendence of the DIS cross section. Hence, the am- 
plitude ^ incorporates geometric scaling. The 
saturation scale Qsix) is parameterized as 



Qs{x) = Qo ( — 
V X 



A/2 



Qo = 1 GeV (2) 



where Xq ~ 3 x 10"^ and A ~ 0.3 [J. 

To investigate whether saturation may be the 
cause of the observed scaling, one can study ge- 



ometric scaling in other experiments where simi- 
larly small values of x are probed, like d-Au scat- 
tering at RHIG. In terms of the dipole amplitude 
the cross section of hadron production in high- 
energy nucleon-nucleus collisions is given by 



K{yh 



dyhd'^Pt (27r) 



Xp 



(3) 



Here Np describes a quark scattering off the 
small- a; field of the nucleus, while Na applies to a 
gluon. The parton distribution functions fq/p and 
the fragmentation functions Dj^/q are taken at the 
scale — pj: which we will always take to be 
larger than 1 GeV^. The momentum fraction of 
the target partons equals X2 = xi exp(— 2yft). To 
good approximation one can neglect finite mass 
effects, i.e. equate the pseu do-rapid ity 77 and the 
rapidity yt and use xp = \Jpi ^ w? /^/s exp(77) « 
Pt / y/s eyip{yh) ■ Finally, there is an overall yh- 
dependent if -factor that effectively accounts for 
NLO corrections. The if-factors are close to 1 in 
the forward region and become relevant towards 
mid-rapidity. 

Unlike in DIS, in hadron-hadron collisions, a 
scaling dipole amplitude will not lead to a scaling 
property of the cross section, due to convolutions 
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with the parton distributions and fragmentation 
functions in ([3]). Hence one cannot establish geo- 
metric scahng in the data directly, but instead one 
has to test the scaling properties of the dipole am- 
plitude using a model, like the GBW model ([T]). 
However, due to its exponential fall-off at large 
transverse momentum, the GBW model cannot 
describe the RHIC data. Instead, a modification 
of the GBW model, which we will refer to as the 
DHJ model, was proposed in f 4|5| . It offers a good 
description of hadron production in the forward 
regioro, and is given by 



NAiqt,X2) 



l-exp(-i(r2Q2(^2))7te..2)) 



(4) 



The corresponding expression Np for quarks is 
obtained from Na by the replacement {r'lQl)'^ — > 
{{Cp/CAylQlV , with Cf/Ca = 4/9. Note that 
the so-called "anomalous dimension" 7 is a func- 
tion of qt rather than rt , so that the Fourier trans- 
form can be obtained more easily. 

The anomalous dimension 7 of the DHJ model 
is parameterized as 



7DHj(gt,a;2) 
+ (l-7s) 



7s 



\\0g{qj/Ql{xM 



\y + d^+\\og{qj/Ql{x2))[ 



(5) 



Here y — log 1 /x2 is minus the rapidity of the 
target partons. The saturation scale Qs{x2) and 
A are taken from the GBW model ([2]), including 
a larger value of Qo ~ 1-63 GeV to account for 
the size of the nucleus. The parameter d — \.2 
was fitted to the data. Away from Q^, 7dhj rises 
towards 1 like {log qt/Qs)/y', clearly the model vi- 
olates geometric scaling. At the saturation scale 
7DHJ assumes a specific value 7^ = 0.628. The 
logarithmic rise, the scaling violations and the 
value of 7s are expected from small- .t evolution 
[B]. However, an analysis of the BK equation 
suggests that a smaller value of 7 w 0.44 may 
be more appropriate at Qs |7]. We note that the 
DHJ model ^ was meant to apply only outside 



-"^As it turns out the study of Ref. contained an artifi- 
cial upper limit on the xi integration to exclude large X2. 
Without this cut, the larger pt data for = 0, 1 are not 
well-described by the DHJ model. 



the saturation region, which is hardly probed at 
RHIC. Hence the behaviour of (O for qt < Qs is 
irrelevant. 

To investigate whether the RHIC data are com- 
patible with geometric scaling, we write down a 
new model that features exact scaling, unlike the 
DHJ model. Further, the new model, which we 
will refer to as the GS model, does not have the 
logarithmic rise or the value of 7 at Qs that are 
both expected from small- a; evolution. Instead, it 
is parameterized as 

7Gs(^=t)=7i + (l-7i) (jr°;)^|^ . (6) 

The parameters a and b will be fitted to the data. 
Not only is 7gs exactly scaling, it also rises much 
faster towards 1 at large qt than 7dhj- Expanding 
the exponential, we see that the large qt behavior 
of Eq.'s dSl) and IH) is given by [8], 
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NA{qt) fx 



.2 + a 
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for 7GS • 



(7) 



Hence, the pt distribution resulting from Eq. ([5]) 
will fall-off much faster than in the DHJ model. 
We note that our model ([6]) is not intended to re- 
place other, theoretically better motivated mod- 
els but is constructed simply to investigate in 
a general way which conclusions can really be 
drawn from the RHIC data in the central and 
forward regions. 

In Fig. [1] we show dNh/ [dytd'^pt) following 
from Eq. ([3]), calculated using the anomalous 
dimension 7gs in combination with the am- 
plitude dU. All Pt distributions of produced 
hadrons measured at RHIC in d-Au collisions 
are well described. At Qs we have chosen for 
7Gs(w = 1) = 71 the same value 7s — 0.628 as in 
the DHJ model, for ease of comparison. We also 
use the same parameterization of Qs{x). We ob- 
tain the best fit of the data for a = 2.82 and 
h = 168. As mentioned, this LO analysis re- 
quires a i^-factor to account for NLO corrections, 
which are expected to become more relevant to- 
wards central rapidity. The p^-independent K 
factors we obtain for y/i = 0, 1, 2.2, 3.2, 4 are given 
by if = 3.4,2.9,2.0,1.6,0.7 for our model, and 
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K = 4.3, 3.3, 2.3, 1.7, 0.7 for the DHJ model. For 
further details of the calculation see i8j. 




p,[GeV] 

Figure 1. Transverse momentum distribution of 
produced hadrons in d-Au collisions as measured 
at RHIC (symbols) for various rapidities yh- To 
make the plot clearer, the data and the curves for 
= 0, 1 and 2.2 are multiplied with arbitrary 
factors, namely 16, 4 and 2, respectively. 



We can conclude that the RHIC data are com- 
patible with geometric scaling at all rapidities. 
Hence, no scaling violations can be claimed to be 
observed at RHIC. Of course, from this analysis 
such violations cannot be ruled out cither. The 
fact that in the forward region where yh = 2.2 — 4 
both the exactly scaling GS model and the scaling 
violating DHJ model describe the data indicates 
that to reach a conclusion on geometric scaling, 
a larger -y/s is needed so that a larger range of 
rapidities and pt can be probed. Further, the log- 
arithmic rise of 7dhj is ruled out in the central 
region = 0, 1. However, the DHJ model breaks 
down only at mid-rapidity and pt ^ 2.5 GeV, 
where the values of X2 become larger than 0.01. 
Even though small- 2; evolution may not be ex- 
pected to be valid anymore at 2: > 0.01, we note 
that in d-Au collisions at RHIC Qs is still larger 
than in DIS at a; = 0.01. 

As mentioned, in the forward region both mod- 
els ([5]) and ([6|) work, although they have very dif- 
ferent properties. To illustrate how sensitive the 
data are to the behavior of 7, Fig.[2]shows various 
7gs(w')'s that describe the available data equally 



well. All are of the form ([6]), with different val- 
ues of 71, a and b. The parameters that define 
the edges of the allowed region are 71 = 0.50, 
a = 2.60 and h = 70.2 for the upper curve, and 
71 = 0.75, a = 3.10 and h — 451 for the lower 
curve. Furthermore, we add lines representing 
7DHJ © for different rapidities. To do so the 
rapidity of the target parton y needs to be ex- 
pressed in terms of w and y^, see [5] for details. 
As said, the region below Q^, where the parame- 
terization of 7DHJ is not smooth, is hardly probed 
at RHIC. Fig. [2] shows that 7 is much more con- 
strained by the data at large w compared to the 
region close to Qs- The reason for this is that 
around the saturation scale rt = l/Qs, the inte- 
grand in the dipole scattering amplitude (|4]) de- 
pends only weakly on 7. In addition, the for- 
ward data, yh = 3.2 and 4, are essentially sensi- 
tive only to 71, since for kinematic reasons they 
probe the region where w is close to 1. There- 
fore, the rise of 7 with w is effectively constrained 
only by the data for yh = 0, 1. From Fig. [2] we 
see that 7dhj falls outside the band of allowed 
7's at large w, but remains inside the band for 
w < 2 — this corresponds to respectively pt ^ 2.5 
in the central region yh = 0, 1 where the DHJ 
curves deviate from the data, and the entire for- 
ward region where both models work. We note 
that a recent study of RHIC data on the nuclear 
suppression factor i?''^" suggests that the DHJ 
model already breaks down at yh = 3.2 [5]. 

Where the DHJ model curves deviate from the 
RHIC data, the probed X2 values are not very 
small and one may argue that a small-x descrip- 
tion cannot be expected to apply in the first place. 
However, at LHC, due to the higher energies, the 
region of small X2 extends to a much larger range 
of Pt , so that the predictions of the DHJ model 
and the new one will be different even at small X2- 
Fig.[3]shows the pt distribution of hadron produc- 
tion in p-Pb scattering at 8.8 TeV at LHC. We 
emphasize again that the prediction of the pt dis- 
tribution using the GS model is to be considered 
a tool for checking whether certain small- a; prop- 
erties are present in the data. For the scaling 
curves the best fit obtained from the RHIC data 
was used. 

For small pt the predictions of the two mod- 
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Figure 2. Various fits of 7Gs(''i'), which describe 
the RHIC data equally well. For comparison we 
show curves representing 7dhj(w, y{w,yh)) at dif- 
ferent rapidities yh- 



els are comparable, since there only the region of 
small values of w is probed where 7dhj and 7gs 
are similar, cf. Fig. O Also in the very forward 
region, i.e. y/j « 6 — 8, only this region is tested 
and one obtains similar results from both models. 

However, there is quite a large range where the 
probed values of X2 are small but the predictions 
are clearly different. From Fig. [3] we see that this 
is the case for rapidities smaller than about 5 and 
Pt > 5 GeV — this region corresponds to approxi- 
mately w >3, where according to Fig.[5]the DHJ 
model indeed deviates from the experimentally al- 
lowed band. Due to the larger energy, the values 
of X2 at LHC are much smaller than in the cor- 
responding region of pt > 2.5 GeV at RHIC. At 
LHC, X2 remains below 0.001 in the entire range 
depicted in Fig. [31 i.e. even in the central region. 
Hence, if the LHC data turns out to be described 
by the GS model, one would have a clear indica- 
tion of consistency with d-Au data from RHIC 
at forward and mid-rapidity within the small-x 
framework of Eq. ^ . 

The main difference between the GS model and 
the DHJ model is the slope of the resulting pt dis- 
tribution, which is directly related to the rise of 
7 towards 1. Hence, a measurement of the slopes 
of the Pt distribution at moderate rapidities yh 
at LHC allows a discrimination between the two 



models in a region where small- a; physics — and 
hence a description in terms of Eq. ([3]) — is ex- 
pected to be applicable. The slower fall-off of the 
Pt distribution in the DHJ model is caused by 
the logarithmic rise of 7 towards 1, which is a 
generic signature of small- a; evolution. Hence, if 
the Pt distribution of the produced hadrons falls 
off much faster than predicted by the DHJ model, 
current expectations of small x evolution do not 
hold at LHC. 
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Figure 3. Hadron production cross section in p- 
Pb collisions at LHC from 7gs and 7dhj at dif- 
ferent rapidities yh- 
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